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Abstract 

We apply linear-response analysis to the Gross-Pitaevskii equation to obtain 
the excitation frequencies of a Bose-Einstein condensate in a vortex state, and 
apply it to a system of Rubidium atoms confined in a time-averaged orbiting 
potential trap. The excitation frequencies of a vortex differ significantly from 
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those of the ground state, and may therefore be used to obtain a spectroscopic 

signature of the presence of a vortex state. 
PACS Numbers: 3.75.Fi, 67.40.Db, 67.90,+z 
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The recent attainment of quantum degeneracy conditions in magnetically-trapped, dilute 
alkali vapors has opened an avenue for studying the many-body physics of Bose-Einstein 
condensates (BECs) in unprecedented detail. Recent experiments have mapped out many 
of the basic properties of alkali BECs: the critical temperature T c the temperature 
dependence of the condensate fraction N /N ||, the contribution of particle interactions to 
the internal energy (i.e. departures from ideal-gas behavior) ||, an d the energies of low-lying 
excitations |5|||. 

There remain many aspects of "classical" superfluid behavior that have not yet been en- 
countered in the atomic BECs. The one we discuss in this paper is the generation of vortices, 
which to our knowledge have not yet been observed in trapped gases. Recent theoretical in- 
vestigations of axially-symmetric harmonic traps have identified vortex states of condensates 
which have sharp values of the azimuthal component m of the angular momentum J7| [| . 
Rotation of the trap at a critical frequency iOcnt, which is of the order of the harmonic trap 
frequency u, should force the condensate into the vortex state [0-0. Detection of this state 
by current imaging techniques is complicated by present magnet geometries, which constrain 
viewing of a condensate to be done more or less perpendicular to the trap symmetry axis 
(hereafter the z axis). Most schemes acquire images of the condensate density integrated 
along the line of sight, and the integrated density of a vortex state perpendicular to the 
trap axis is not much different from that of the condensate ground state. In this paper 
we calculate the mechanical excitation spectrum of a vortex state, and find that it differs 
significantly from that of the ground state. We thus propose excitation spectroscopy as a 
sensitive technique for detecting the presence of vortices. 

The basic framework of our method is mean-field theory as described by the Gross- 



Pitaevski (GP) formalism [|I0| , |TT| for a condensate of a dilute Bose gas at temperature T = 0. 
Current experimental BEC atomic physics appears to be practiced in a regime where this 
formalism is applicable: the gases are very tenuous, and a nearly pure condensate (cor- 
responding to T very near zero) can be produced by the technique of evaporative cooling 
]l]Jl]|l]. Calculations done within the GP framework yield good agreement with experiment 



concerning condensate shapes and sizes 0^,0,0, and give ground-state condensate exci- 
tation frequencies within about 5% of experimental values PJ6|fl^]. 

We begin with the time-independent treatment of condensate eigenfunctions using the 
Gross- P it aevskii (nonlinear Schrodinger) equation, and then calculate the excitation spec- 
trum using the method of Bogoliubov [15], which was used by Pitaevskii to examine exci- 



tations about vortices in a homogeneous gas fl9|] , and has been adapted to treat trapped 



Bose condensates [|l4l , |l6Hlq1 ■ The specific calculations reported here are done for a system 
of 87 Rb atoms confined in the time-averaged orbiting potential (TOP) trap currently in use 



at JILA HQ 



In the GP formalism the interaction between atoms is approximated by the pseudopo- 
tential, V(r, r') = U 5(r - r'), where U = 4irh 2 a/M, with M being the atomic mass and a 
the scattering length that characterizes low-energy atomic collisions. For the triplet Born- 
Oppenheimer potential of 87 Rb2 that describes collisions in the JILA trap, the current best 
estimate JH] of a is 110a , where a is the Bohr radius. This value is used in the present 
paper. The trapping potential takes the form V^ rap (r) = M(u 2 p 2 + u 2 z 2 ) /2, in the cylindrical 
coordinates appropriate to the TOP trap. For this trap, the axial and radial frequencies are 
always in the same ratio, uj z /uj p = \/8; it is convenient to characterize the trap by the single 
parameter v p = uj p /2-k, the value of which is often stated explicitly in experimental papers. 
In this paper we use v p = 74 Hz. 

The time-independent GP equation that describes the condensate wavefunction ip(r) 
thus takes the form 

/ h 2 \ 

V 2 + K rap (r) + N U \^(r)\ 2 V(r) = /^(r), (1) 



V 2M 

where N is the average number of condensate atoms, and \x is the chemical potential, which 
is treated as an eigenvalue; the normalization condition 



dr|^(r)| 2 = l (2) 
is implied. The form of Eq. (|l]) is consistent with the existence of solutions 



i>(r)=^(p,z)-= (3) 



that are eigenfunctions (with eigenvalues hm) of the azimuthal angular momentum operator 
l z . Previous work on this system |T|^ |l^|l~3f has identified a solution with m = as the 
condensate ground state, and those with \m\ = 1 as the lowest vortex solutions. 

As discussed elsewhere P^fl , we solve Eq. (|T|) by "growing" a condensate up from the 
noninteracting case (Uq = 0), in which the solution is given by harmonic oscillator wave 
functions. The growth process may be visualized as a gradual turning-on of the interaction 
strength Uq or as a gradual increase in the number of condensate atoms N ; the two pictures 
are equivalent for this purpose since N and Uq appear in the equations of motion only 
through the product NqUq. In the noninteracting case, m is a good quantum number, and 
it is preserved during the growth process. We start with tp^ m \p,z) as an eigenfunction of 
the axially symmetric three-dimensional harmonic oscillator with quantum numbers m, the 
number n z of nodes in z, and the number n p of nodes in the cylindrical radial coordinate p. 
The lowest vortex state with m = 1 has n z = n p = 0; the number of nodes in p and z is also 
found to be conserved during the growth process. Our general representation of ip( m \p, z) 
is as a linear combination of oscillator eigenfunctions, and a set of nonlinear equations is 
solved iteratively to determine the coefficients at each value of N U . A cross-section of the 
density of the vortex solution is shown in Fig. |l|. 

With the solution of the time-independent GP equation in hand, we calculate the re- 
sponse of the condensate to a weak oscillatory perturbation by standard linear-response 
theory p3 |. The associated time-dependent driven GP equation takes the form: 



ih 'dt 



- #o + WM)| 2 

+ f(r)e- luJpt + /*(r)e^ t l (r, t), (4) 



where /(r) is the spatially-dependent amplitude of the perturbation. We solve this equation 
in the linear- response limit. The details of this approach are described elsewhere |TBJ, and 
we simply state the central results here. By using the form 



*(r,t) = e 



-na/n N^{r)+u{r)e- itJpt + v*(r)e iuj » t (5) 
we obtain Eq. ([]]) and also the linear-response equations, 

(£ - fiu; p ) u(r) + N U mr)] 2 v(r) = -iV 1/2 /(r)^(r), (6) 

N U ir(r)} 2 u(r) + (£ + ^>(r) = -iV 1/2 /(r)^*(r), (7) 

where £ = H - fx + 2N U |^(r)| 2 . 

This pair of equations can be solved in a general way by expanding the condensate 
response in normal modes of oscillation, which are obtained by solving the Bogoliubov 
equations, 

(£ - hu x ) u x (r) + iVWo [^(r)] 2 v x (r) = 0, (8) 

Wo [r(r)] 2 u x (r) + (£ + fo; A ) v A (r) = 0, (9) 

where u; A is an eigenvalue and u x (r), V\(r) are corresponding eigenfunctions. We require that 
u\(r) and v\(r) be square integrable and satisfy the conventional normalization condition 

J dr (\u x (r)\ 2 - \v x (r)\ 2 ) = 1. (10) 

With this condition in force, the condensate response to an arbitrary driver /(r) is obtained 
by a superposition of normal modes , 



W =-y H r) ^) (in 

where the amplitudes g x are obtained by quadrature: 

g x = N 1 /2 J dvf(v)[^(v)u* x (r)+r(r)v* x (r)]. (12) 

Two key qualitative aspects of condensate excitation are implied by Eqs. ( |TTD and (P"2"|). 
First, from Eq. (|TTJ) , the response is largest when the driving frequency u p is equal to a 
normal-mode frequency uj x (the apparent divergence in this response, as in the case of a 



driven classical oscillator, is due to the neglect of damping in this treatment). Second, Eq. 
(|T2|) implies selection rules for a given driver /(r), associated with symmetries of the solutions 
of Eqs. (H) and (|9]). In particular, if ip{r) is given by Eq. (|3|), then it is straightforward 
to show that the normal modes will have specific angular momentum composition, in the 
following sense: if U\(r) is an eigenfunction of l z with a particular eigenvalue m u (in units 
of h), then v\(r) will be an eigenfunction with eigenvalue m u — 2m. It is appropriate to 
think of a normal mode as a quasiparticle moving in an effective potential created by the 
condensate, and for v\ as being created by scattering of u\ by the condensate. For m — 0, 
the condensate has axial symmetry, and the normal modes thus have definite values m u 
of the angular momentum; for m 7^ 0, the condensate is not axially symmetric, and the 
component U\ (with angular momentum m u ) is scattered into v\ (with angular momentum 
m u — 2m). In the case treated in this paper, we have m = 1, and we will label the normal 
modes by the value of m u that corresponds to the component U\. Thus, a "breathing-mode" 
driver, such as was applied to the ground state of the condensate in Ref . |j| , will result here 
in excitation of dipole driver, upon which we comment below, yields m u = and 

2; and a quadrupole driver, also used in Ref. |5[], gives m u = — 1 and 3. 

We have solved Eqs. (|] ) and (^) by an extension of the technique used in Ref. JTJ 
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which u\ and V\ are expanded in trap eigenfunctions of appropriate symmetry, to obtain a 
system of linear eigenvalue equations. Fig. [I] depicts ip along with u\ and V\ as computed 
for the lowest normal mode, which has m u = 2, and would therefore be generated from 
the vortex condensate by a dipole excitation (e.g. oscillatory displacement of the center of 
the trap). The s- and <i-wave characteristics of v\ and u\ respectively are apparent in the 
figure. An important property of this normal mode, which differs from all cases we have 
encountered in excitation of m — 0, is that its frequency uj\ is less than the trap frequency 
Up. The eigenfrequency of the lowest dipole mode for m = 0, on the other hand, is identical 
to the trap frequency; that mode describes the rigid oscillation of the condensate ground 
state. 

The dependence of normal mode frequencies upon iVo is depicted in Fig. |^. It shows 



two candidates for excitation by dipole or quadrupole driving that have distinctive spectral 
signatures, in that no excitation frequencies of the m = condensate are nearby. In addition 
to the m u = 2 mode just described, there is the m u = —1 mode, which could be excited by 
a quadrupole rotation in the sense opposite to that of the vortex. This frequency vanishes 
in the noninteracting limit, due to the m degeneracy of the cylindrical harmonic oscillator; 
its non-zero value elsewhere results directly from interactions of the quasiparticle with the 
condensate. 

It is worth noting that the frequency spectrum in Fig. |] is applicable to a wide range of 
TOP-trap geometries. Solutions of Eqs. (§) and @, subject to the constraint of Eq. (|10|), 
apply to all TOP-trap geometries for which the parameter 

1 = N a(Mu p ) 1/2 (13) 



remains invariant |fL4|j . This scaling law makes excitation spectroscopy a particularly effective 
tool for diagnosing the presence of vortex condensates. 

In conclusion, we have calculated the excitation spectrum of vortex states of dilute Bose- 
Einstein condensates, in the linear response regime. These spectra are sufficiently different 
from those for ground state excitation, that they may provide useful diagnostics of the 
presence of vortices in trapped atom systems. 
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FIGURES 

FIG. 1. The dashed line is a plot of 100 \ v\ (r)| 2 for the lowest energy excitation, the dotted 
line shows 100 \u\ (r)| 2 , the factor of 100 is an arbitrary scaling factor used for plotting convenience. 
The solid line shows a plot of the spatial distribution of the k = 1 vortex number density, n (r). 
In this plot N = 4939 atoms, and v p = 74 Hz for a TOP trap. 

FIG. 2. Normal-mode excitation frequencies oo\ of the vortex state in the TOP trap vs. number 
of 8 Rb condensate atoms No, in units of the trap frequency w p . Labels indicate the angular mo- 
mentum quantum number m u of the u\ component of the normal-mode eigenfunction, as described 
in text. There are two degenerate frequencies at lo\ = oj p for all values of iVo, which correspond to 
rigid oscillations of the vortex in the trap. 
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R. J. Dodd et al. Fig. 2 



